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I . INTRODUCTION 


I 2 

In  the  theory  of  the  elastic-perfectly  plastic  solid  * two 
distinct  types  of  constitutive  equations  are  usually  applied  to  the 
same  material:  one  type  is  used  in  the  elastic  range  and  the  other 
type  is  used  in  the  plastic  range.  A similar  statement  could  be  made 
for  an  elastic-plastic  strain-hardening  material.  Our  purpose  here  is 
to  show,  by  example,  that  a single  set  of  equations  of  the  rate  type 
can  describe  behavior  both  in  the  elastic  and  in  the  plastic  range. 
Furthermore,  these  equations  can  describe  both  loading  and  unloading. 
Indeed,  the  mathematical  machinery  needed  to  describe  elastic-plastic 
behavior  is  contained  in  the  classical  theory  of  ordinary  differential 
equations.  We  shall  show  how  elastic  behavior  can  correspond  to 
uniqueness  of  solutions  of  such  equations;  how  nonuniqueness  of  solution 
can  correspond  to  plastic  flow,  and  how  the  possibility  of  different 
constants  of  integration  can  explain  why  an  elastic  unloading  curve  can 
differ  from  a previous  elastic  loading  curve. 

II.  A ONE-DIMENSIONAL  EXAMPLE 

Before  proceeding  to  our  proper  three  dimensional  example,  we  shall 
present  a one-dimensional  example  in  order  to  exhibit  how  some  of  the 
mathematical  notions  from  ordinary  differential  equations  theory  apply. 

We  consider  a single,  strain  component  y and  a single  non-dimensionalized 
stress  component  s which  will  correspond  in  our  three-dimensional  example 
to  shear  strain  and  shear  stress  respectively  in  a simple  shearing  motion. 
We  assume  that  stress-rate  is  related  to  strain-rate  and  stress  through 
an  equation  of  the  form 


ds  _ 1^  L _n  2n  ^ 

dt  “ 2 ^ s dt  ’ 

where  t is  time  and  n is  some  positive  integer.* 


(1) 


have  chosen  to  work  in  terms  of  non-dimensional  stress  in  order  to 
avoid  the  necessity  of  defining  symbols  for  material  constants^  which 
the  reader  then  is  required  to  keep  straight.  A reformulation  in  terms 
of  such  constants  can  readily  he  made  when  needed. 

^Prager,  W.  and  Hodgej  P.  G. , "Theory  of  Perfectly  Plastic  Solidss" 
Dover,  New  York,  L.  C.  68-19164  (1968). 
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Thomas,  T.  T.,  "Plastic  Flow  and  Fracture  in  Solids,"  Academic  Press, 
New  York,  L.  C.  if  61-12277  (1961). 


Preceding  Page  Blank 


In  the  one-dimensional  case  currently  under  discussiont  we  shall  find 
it  convenient  to  eliminate  t and  write  merely  for  (1) 


37  ■ T . C2) 

For  simplicity,  we  shall  first  discuss  the  case  n = 1.  In  this  case  the 
solutions  of  (2)  are 


The  restrictions  on  the  range  of  y in  which  the  solutions  (3)  are 
valid  are  necessary  since  (2)  requires  that  ds/dy  be  non-negative.*  The 
solutions  (4),  in  which  s takes  one  of  two  constant  values,  may  seem 
rather  special.  However,  they  play  an  essential  role  in  our  formulation. 

Imagine  now  that  starting  with  y <=  0,  s > 0,  one  begins  to  increase 
y.  The  only  solution  of  (4)  (with  n - 1)  possible  with  these  initial 
conditions  is 


s « _i  sin  -I  . (S) 

ii 

As  Y is  increased,  s will  be  given  by  (5)  and  will  continue  to 
be  given  by  (5)  as  long  as  y remains  less  than  'll  ti/2.  Indeed,  the 
Piccard-Lindelof  uniqueness  theorem^  assures  us  of  this,  since  a 
Lipschitz  condition  will  hold  when  -1//2  < s < 1//2.  Indeed,  as  long 
as  |s|  does  not  reach  1//2  , s will  follow  the  elastic  solution  (S) 
whether  y is  increased,  decreased  or  otherwise  varied. 


Of  Qourse  we  aan  just  as  well  let  y//S  lie  in  a Tange 
2m  + a - ~ " — 2m  + cj  + j foT  some  integer  m other  than  zero. 

V 2 

However  this  would  just  introduce  an  unneaesaary  redundancy,  since  Hum 
aan  be  abaorbed  into  the  constant  a. 

^Coddingion,  E,  /}.  and  Levinson,  N.,  ^'Theory  of  Ordinary  Differential 
Equationa,"  McGraw-Hill,  New  York  (1966). 
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Now  imagine  that  y is  increased  until  it  reaches  the  value  in/2/2 
and  is  thenceforth  increased  further.  The  stress  will  have  then  reached 
the  value  1//2  at  which  the  Lipschitz  condition  fails,  and  so  uniqueness 
will  no  longer  be  implied  by  the  Piccard-Lindelof  theorem.  Since  no 
real  solutions  of  (2)  (with  n » 1)  exist  for  s > 1//2,  and  since  s cannot 
decrease  when  y fs  increasing,  it  follows  that  s will  remain  at  the  value 
1/i/T  as  long  as  y is  not  decreasing.  In  other  words,  in  loading,  after 
s has  reached  the  first  solution  (4)  takes  over.  This  is  a plastic 

yield  solution,  and  it  must  prevail  at  the  critical  value  of  s,  namely 
s = 1//2.  This  process  is  illustrated  in  Figure  1. 

Now  suppose  that  yielding  proceeds  as  described  above  until  the 
strain  has  reached  the  value  y,  > •n/Ell.  Suppose  that  at  this  point 
the  strain  begins  to  decrease.  The  differential  equation  (2)  then  allows 
the  stress  to  follow  more  than  one  solution:  It  allows  the  solution 


^ y-yj^  + ti/Z/l 

s = — sin  


(6) 


but  it  now  also  allows  the  solution 


s = 1/*^  . (7) 

since  both  for  (6)  and  (7) , ds/dy  ^ 0 when  y is  decreasing  and  s is 
not  increasing.  However,  note  that  for  decreasing  y,  (7)  is  an  asympto- 
tically unstable  solution  in  the  sense  that  a small  perturbation  in  s will 
cause  the  solution  to  move  onto  one  of  the  forms  (3)  which  then  takes  s 
far  from  1//2  as  y decreases  further.  However,  (6)  is  neutrally  stable 
in  the  sense  that  a small  perturbation  in  s will  not  cause  an  increasing 
divergence  of  the  perturbed  solution  from  (6).  If  we  make  the  assumption, 
then,  that  the  stable  solution  will  be  chosen,  then  we  see  that  s will 
follow  the  solution  (6),  as  shown  in  Figure  2.  As  soon  as  s is  below 
1/*^,  the  (s,y)  values  will  have  again  entered  the  region  of  uniqueness 
of  the  differential  equation  (2)  and,  thus,  will  stay  on  the  curve  (6) 
whether  y is  decreased  or  even  increased  again  as  long  as  the  increase 
or  decrease  is  not  so  largo  as  to  bring  s to  one  of  the  yield  values  y 
= ± 1/*^.  There  is  a return  to  elastic  behavior  about  a new  stress  free 
state,  that  for  which  y = y^  - ■(t/2. 

In  OUT  formulation  then,  elaetio  behavior  ia  thua  related  to  unique- 
neaa  of  aolutiona  of  a differential  equation.  Different  elaatia  regimes 
aorreapond  to  different  aonatanta  of  integration,  Plaatia  behavior  ia 
related  to  nonuniqueneaa. 

Let  us  now  go  a bit  further  with  our  example.  Suppose  that  we 
continue  decreasing  till  we  reach  s •=  - Again  there  will  be  plastic 

yield  (see  Figure  3) ; s cannot  decrease  further  and  cannot  increase  as 
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long  as  Y is  decreasing.  Thus  s will  follow  The  solution  s = -\l /Z 
until  Y is  made  to  increase  again,  after  which  it  will  follow  another 
one  of  the  solutions  (3).  It  will  enter  another  elastic  regime.  If  y 
is  increased  and  then  decreased  in  a periodic  manner  over  several  periods 
beyond  the  point  that  s reaches  the  critical  value  \//2  on  each  increase 
and  -\l</2  on  each  decrease,  the  solutions  will  follow  a hysteresis  loop 
(see  Figure  4) . 

Now  one  objection  to  the  above  argument  is  that  in  the  elastic 
regimes  the  (y,  s)  diagram  follows  a sine  curve,  whereas  it  would  be 
preferable  to  have  it  follow  a straight  line  until  it  is  near  yield  and 
then  curve  abruptly  to  join  the  yield  curve,  say  s » 1/v^.  Such  behavior 
can  be  realized  by  choosing  a large  enough  value  of  n,  as  we  shall 
illustrate  now.  Indeed,  solutions  to  (2)  are  given  by 


(8) 


where  B^(a,b)  is  the  incomplete  beta  function,  and  where  the  plus  sign 
is  used^if  s and  the  minus  sign  if  s <^0.  (The  function  so  defined 
is  analytic  for  -1  < s < 1).  Yield  occurs  at  s=>±l//2> 


Indeed,  the  behavior  is  qualitatively  the  same  as  in  the  example 
for  n = 1,  except  that  linear  elastic-perfectly  plastic  behavior  joined 
by  a transition  knee  is  approximated  for  largo  n. 

In  Figures  S,  6 and  7 we  show  curves  of  stress  versus  strain  increas 
ing  beyond  yield  for  n = 16,  32  and  64. 


III.  A THREE-DIMENSIONAL  EXAMPLE  - ELASTIC-PERFECTLY  PLASTIC 
SOLID-PRELIMINARY  ANALYSIS 


We  shall  treat  here  a three-dimensional  example  which  corresponds 
to  an  elastic-perfectly  plastic  solid  undergoing  infinitesmal  strain. 

The  solid  will  have  a strain  energy  function  valid  for  each  of  its 
elastic  regimes.  It  will  exhibit  yield  in  the  form  of  nonuniqueness 
of  solutions  of  its  governing  rate-type  constitutive  equations.  And, 
again,  elastic  behavior  will  be  associated  with  uniqueness  of  solutions 
of  ordinary  differential  equations  and  plastic  behavior  with  nonuniqueness . 
Moreover,  in  the  plastic  regime,  it  will  follow  the  Prandtl-Reuss  equations, 
not  as  a new  assumption,  but  as  a consequence  of  the  governing  differential 
equations . 


4 

AbrojnowitZj  M,  and  Stegunj  I.  A.,  "Handbook  of  Mathematiaal  Funationst" 
U.b.  Government  Printing  Office,  Washington,  D.C.,  L.  C.  if  64-600S6 
(1964). 
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Thus  (15)  and  (16)  yield 


■ 3 '^kS.'^ijj 


da 

= “T? 


We  obtain,  then,  from  (13),  (14),  and  (17), 


a.  . = Ke,  • ■*■  £•  • . 

kk  ij  a 13  ’ 


whence,  using  (9)  and  (10),  we  have 


°kk  = ^‘^^kk 


^ ... 

IJ  a 13 


Squaring  both  sides  of  (20),  summing  and  applying  (12)  and  (13), 
wo  got 


6^=4..4..»  e.  . e.  . = [(i»’(a)1^ 

13  13  ^2  13  13  i.y  V ; j 


6 = (!>'  (a)  . 


(Note,  wo  have  assumed  that  i}i'(a)  > 0 in  some  neighborhood  of  a = 0 so 
that  (22)  will  follow  from  (21)  there.  Remember  that  B > 0) . 

Lot  the  solution  to  (22)  be  given  by 


a = i|)  3 


.1 


Then  from  [20) , (22) , and  (23) , we  get 


•4. . 


ij  'i'CP)  ij 


(24) 


or 


i^s..  = e.. 

6 13  13 


(25) 


Differentiation  of  (25)  with  respect  to  time  (where  dot  denotes  time 
differentiation)  yields 


4 ^ m 4 = *£. 

e i3  fl2  ij  s ij  Ij 


(26) 


Multiply  both  sides  of  (26)  by  4^^  and  sum  to  get 


Now  (12)  yields  by  differentiation 


0B  = .4. , 

13  13 


(28) 


Use  (12)  and  (28)  in  (27)  to  obtain 


or 


0 = 


0i|)'(B) 


(29) 
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If  we  now  substitute  the  right  hand  side  of  (29)  for  0 in  (26)  and 


then  solve  for 


we  obtain 


^ij  " ij;(e) 


I 


(30) 


At  this  juncture,  we  make  a fundamental  new  assumption  which  is  the 
basis  of  our  departure.  Vie  assume  that  the  basic  constitutive  equations 
are  (19)  and  (30)  instead  of  (19)  and  (20).  We  oonsidev  a material  of  the 
rate  type  for  which  the  differential  equations  (30)  are  the  governing 
aonetiUitive  relations  between  stress  deviator  and  strain  deviator. 


Thus  behavior  described  by  (20)  can  occur  as  one  possible  solution 
of  (30).  Other  solutions  are  obtained  by  replacing  by  + Cj^^ 

in  (20),  where  Cj^^  is  a traceless  constant  tensor:  this  can  be  seen 

immediately  from  (30)  since  the  time  derivatives  of  and 

are  the  same  and  the  strain  deviator  enters  (30)  only  through  its 
derivatives. 


Furthermore,  equation  (29),  which  can  be  derived  from  (30)  by 
multiplying  both  sides  by  and  summing,  plays  a critical  role  in  the 

transition  to  the  plastic  regime.  We  shall  now  show,  by  specific  examples, 
how  with  certain  choices  of  function  i|>(6),  the  equation  (29)  will  govern 
the  onset  of  plastic  yield. 

IV.  THREE- DIMENvSIONAL  EXAMPLES  OF  ELASTIC -PERFECTLY  PLASTIC 
SOLIDS  GOVERNED  BY  A SINGLE  EQUATION 

We  now  seek  examples  of  functions  iji,  such  that  equation  (30)  will 
give  olastic-perfectly  plastic  behavior.  Please  note  that  although  in 
our  derivation  in  Section  III  above,  we  used  (29)  to  get  to  (30),  in 
fact  it  is  not  difficult  to  establish  that  (29)  also  follows  from  (30). 

Now  let  us  observe  that  if  we  had 


(31) 


14 


then  (29)  would  read 


(32) 


2 

This  would  then  imply  the  von  Mises  yield  condition.  Indeed  if  g 
rises  to  the  value  unity,  then  (32)  implies  that  as  long  as  positive 
work  is  being  done  by  the  shearing  stresses,  cannot  decrease  and 
thus  must  remain  at  unity. 

Now  (31)  gives  us  a differential  equation  for  >(»(&).  Integration  of 
this  equation,  with  the  constant  of  integration  chosen  so  that  i(/(g)/B 
has  a removable  singularity,  and  is  thus  well  defined,  at  g “ 0 gives 


)Jj(g)  - sin"  g 


(33) 


(To  see  why  we  want  iii(g)/g  to  be  well  defined  at  6 = 0,  see  how 
this  expression  figures  in  (25)) . 

If  (33)  is  then  used  for  ij),  (30)  becomes 


sin'^g 


11 


(34) 


Observe  then,  that  for  an  assignment  of  as  a continuously 

differentiable  function  of  time  (or  even  as  a piecewise  smooth  function 
of  time)  the  Piccard-Lindelof  theorem  implies  a unique  solution  3..(t) 


,2 


ij 


of  (34)  provided  g <1.  If  as  remarked  above,  g reaches  unity  with 
positive  shear  work  continuing  to  be  put  in,  then  g^  remains  equal  to 
unity,  and  we  have  plastic  flow  with  a von  Mises  yield  condition.  During 
plastic,  yield,  we  put  g = 1 in  (34)  to  obtain 


ij 


IT 


4- .4,  „ 
kX. 


(35) 


Note  that  the  equations  (35)  are  the  Prandtt-Reuse  equations. 
(The  reader  is  reminded  that  the  material  constants  have  been  given 
specific  values  for  simplicity.  Otherwise  they  should  appear  with 
arbitrary  values  in  the  Prandtl-Reuss  equations) . 
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Next  let  us  look  at  the  solutions  to  (34)  in  the  region  of 
uniqueness,  < i,  (i.e..  Let  us  look  at  solutions  in  the  elastic 
regime.)  If  we  compare  (33)  and  (23)  and  (22),  we  get 


3 = sin  0=4'  (“) 


(36) 


whence 


4(01)  = 1 - cos  a,  (37) 

if  we  choose  the  constant  of  integration  of  the  equation  (36)  for  4 
so  that  4(0)  = 0.  (Since  4(a)  figures  in  the  strain  energy  function, 
to  which  one  can  add  an  arbitrary  constant  without  affecting  the  stress - 
strain  relations,  this  choice  of  the  constant  of  integration  is  merely 
a matter  of  taste:  it  insures  that  4(“)  gives  a contiibution  to  the 
strain- energy  which  is  non  negative  and  vanishes  at  a = 0.) 

Now  (13)  with  the  aid  of  (37)  gives 

W = I f^®kk^^  + 1 - cos  a (38) 

whence  (20)  yields 


4 


ij 


sin  a 
a *ij 


(39) 


Now,  (39)  is  a particular  solution  of  (34) . To  get  the  general 
elastic  solution,  we  can,  according  to  the  remarks  at  the  end  of  Section 
III,  replace  by  + c^^  where  c^^^  is  a traceless  constant  tensor. 

Thus,  we  obtain  for  the  elastic  solutions  of  (34) 


4,  . = (e.  . + c.  .) 


(40) 


where  Cj^j^  = 0.  Furthermore,  with  the  restriction  3 ^1,  (36)  yields 
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as  a constraint  on  the  region  of  validity  of  (40)  for  given 
strain  energy  corresponding  to  a given  elastic  regime  is  then 


The 


“ I - 1 - cos  *'C  ^ + Cy.^)  . (41) 

Thus  for  each  elastic  regime,  we  have  a strain  energy  function  (41). 

The  onset  of  plasticity  is  governed  by  equation  (32) . During  plastic 
flow  we  have  the  Prandtl-Reuss  equations  (35).  And  all  of  these  conclu- 
sions are  consequences  of  the  single  set  of  constitutive  equations  (34) . 

We  now  consider  simple  shear  for  which  4^2  “ ^2^  ■ 4,  4^^^  ^ 0 
otherwise,  *12  ” *21  ” 2^  * ^ otherwise.  We  obtain  from  (12) 

02  . ,2^  . ,2^  = 242.  (42) 


We  also  find  that 


4 


ij 


(43) 


If  we  put  (42)  and  (43)  into  (32) , we  obtain  equation  (1)  with 
n 's  1,  and  thus  the  example  of  Section  11. 

For  other  values  of  n,  let  us  put 


♦ '(H  - . ‘ 

^1- 

into  (30) . We  then  get  from 

♦«)  ■ “ (k  ■ i)  • '«) 


For  this  our  analysis  holds  in  a manner  quite  similar  to  that 
developed  above.  Equation  (29)  gives 


d 

dt 


2 


(46) 
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so  that  again  3 cannot  increase  beyond  unity  as  long  as  positive  shear 
work  is  being  done  and  a von  Mises  yield  condition  holds.  During  yield 
one  obtains  from  (30)  and  (44)  the  Prandti-Reuss  equations 


■^ij  ” ij)(l)  ^ij  ' i|)(l)  '*ij  \l  %.!L  * 

During  an  elastic  regime,  one  gets  a stress-strain  relation  of  the  form 


ii ‘‘ij  * “ij’  ■ 

"la  ^ki  'kt 

where  Cj^j,  = 0 and  is  restricted  to  the  domain 


(48) 


and  (j)(a)  is  obtained  from  t;/(3)  by  the  relations  (22)  and  (23).  The  stress 
versus  strain  in  a simple  shear  regime  is  shown  in  Figures  5,  6,  and  7 
for  different  values  of  n.  The  inverse  of  the  elastic  modulus,  namely 
\Ji(3)/3  is  shown  plotted  versus  3 for  various  values  of  n in  Figures  8, 

9 and  10.  Note  that  for  n = 16,  our  calculations  show  that  to  five 
decimal  place  accuracy,  it  does  not  deviate  from  unity  until  6 exceeds 
0.8  and  it  achieves  a maximum  value  of  1.0411  at  yield,  i.e.,  at  6 ® 1. 

(It  is  constant  to  within  4.2%  through  the  whole  elastic  range.)  For 
n = 32,  the  inverse  of  the  modulus  remains  one  to  five  decimal  places 
till  3 > .92  and  it  has  the  value  1.0207  at  yield.  For  n •=  64  it  remains 
unity  to  five  decimal  places  until  3 > .96,  and  it  has  a yield  value  of 
1.0104.  (It  is  constant  to  within  about  1%  in  the  whole  elastic  range.) 

We  see  then  that  with  large  enough  values  of  n,  one  may  approximate  as 
closely  as  one  wishes  the  type  of  elastic-perfectly  plastic  behavior  in 
which  there  is  a sudden  break  in  the  slope  of  the  stress-strain  curve 
at  yield. 


V.  STRAIN  IIARDENING,  AN  EXAMPLE 

We  shall  now  discuss  the  extension  of  the  foregoing  ideas  to  a 
strain-hardening  material.  To  this  end,  consider  any  one  of  the  functions 
<j)(a)  discussed  above  apd  replace  (13)  by 


K.  .2  k2 
" I ^®kk^  + ^ + <}’(«) 


(50) 


where  k is  a material  constant. 
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(For  exaople,  we  may  take  ^(a)  -I  - ccs  a in  (50).  Then,  using  (14), 
instead  of  (20) , we  obtain 


4 = Ke.  + e 

ij  a 13 


(51) 


If  we  define  by 


"ij  ■ •‘ij  - “ij  • 


(52) 


then  we  have  from  (51)  and  (52) 


p . = e.  , 

^3  a 13 


(53) 


If  we  now  define  c by 


C ■ p.  , B,  . 
^13  ^i3 


(54) 


we  obtain  instead  of  (30)  the  equation 


ij  ’j'CC)  ij 


2 [ij)'  * I"]  ^ij  ^ki,  ‘kk 


(55) 


or 


■^ij  “ "^^ij  t ^ij  ^ 2 fijf'  ~ i|)l  ^^ij  * '^^ij^^^kk  " '^'kk^  ' 

^ t rrc. 


We  obtain  instead  of  (29) 

* _ 1 . 

^ cFTcT  ^kk  ‘kk  ' 

If  i))(c)  is  taken  such  that 


'kk- 
(56) 


(57) 


»l''(0  - (1  - 


(58) 
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2 

we  see  that  we  must  have  c and  yield  occurs  when  ? ± 1 . During 

yield  we  get  the  flow  rule 


^ij  " * Ijltrr  ■ itiy  ■ •' W ‘kf 


Uniqueness  of  solutions  of  (56)  holds  for  < 1.  These  solutions 
are  of  the  form 


4 . . a K e.  . 
IJ  1} 


i] 

■^  ( * c^j)  * (60) 


where  Cj^j^  ■ 0.  The  strain  energy  has  the  form 


" ■ J t“kk)^  * y * ♦ (^(“ki  * ‘kiX'ki  ‘ “kij)- 

The  yield  condition  - 1 as  implied  by  (57)  and  (58)  is  then 


C^ij  - Kt..)  (4^,  - KC.p  a 1 


(62) 


For  an  example  of  how  this  differs  from  the  perfectly  plastic  solid, 
lot  us  consider  simple  shear  with  k «=  1/2  and  n = 1.  We  then  obtain 


C » /2  U -1  1 . 

so  that  the  plastic  yield  solutions  are 


The  elastic  solutions  are 


(63) 


(64) 


(65) 
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A plot  of  the  stress  response  to  a strain  in  which  y is  first 
increased  beyond  yield  and  then  decreased  until  yield  again  occurs  is 
shown  in  Figure  11. 

VI.  DISCUSSION  AND  CONCLUSIONS 

Our  purpose  has  been  to  show  by  example  that  elastic-plastic 
behavior  can  be  embodied  in  a single  set  of  constitutive  relations  of 
the  rate  type.  Elasticity  is  then  associated  with  uniqueness  and 
plasticity  with  nonuniqueness  of  solutions  of  the  differential  equations 
which  the  constitutive  relations  comprise.  We  wish  now  to  make  explicit 
a point  which  has  up  to  now  only  been  implicity  stated:  the  plastic 
solutions  do  not  arise  as  asymptotic  solutions  of  the  equations.  They 
are  reached  at  definite  finite  values  of  strain  during  any  given  loading. 
There  is  no  need  to  replace  any  solution  by  an  asymptotic  solution:  the 
plastic  and  the  elastic  solutions  are  exact,  both  for  loading  and  un- 
loading. 

In  closing,  we  ask  the  reader  to  bear  in  mind  that  in  this  report 
we  arc  merely  presenting  examples  to  show  that  it  is  possible  to  unify 
elasticity  and  plasticity  in  a single  rate-type  theory.  Wo  hope  that  in 
the  future  we  can  state  the  notions  herein  contained  more  generally  and 
abstractly  so  that  their  application  can  be  extended  to  more  types  of 
materials  than  those  discussed  herein. 
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Theoretical  Stress-Strain  Curve  During  Loading  Obtained  as  an  Exact  Solution  of  (1) 
with  n = 64  Showing  Elastic  Response  Followed  by  Yield.  Shear  Stress  is  & and  Shear 


Figure  8*  Inverse  of  Half  the  Elastic  Shear  Modulus  as  a Function  of  the  Stress  Invariant  P 

for  ifi  Given  By  (45)  with  n = 16.  Note,it  is  Constant  to  Kithin  4.2%.  Its  Value  at 
Yield  CP  = 1)  being  1.0411. 


Figure  10.  Inverse  of  Half  the  Elastic  Shear  Modulus  as  a Function  of  the  Stress  Invariant  g for 
ii  Given  by  (45)  with  n = 64.  Note  that  it  is  Constant  to  Within  1,1%.  Its  Value  at 
Yield  C3  = 1)  is  1.0104. 
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